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ABSTRACT 


Reactions between atoms and ions of the interstellar 
medium and vibrationally excited molecular hydrogen are 
shown to produce a specific set of molecules. Molecules 
such as CO and CN can readily be produced in abundances 
greater than those of corresponding hydrides. When quantified 
in a simple manner the equilibrium density of molecules 
compares well with the observed values for C Oph. 



I. INTRODUCTION 


We have recently discussed the role played by internal 
energy in the ion-molecule reaction C + 4- Hg "* products, and 
the relevance of this reaction to the interstellar medium 
(Stecher and Williams 1972a). We have argued on the basis 
of many similar examples that sufficient molecular vibrational 
energy can overcome endothermicity and that such reactions 
then proceed with high rate. Our purpose here is to 
discuss in a similar way the chemical exchange reactions 
(sometimes called metathetical reactions). We are then 
able to list a set of many reactions which will be expected 
to proceed with high rates in the interstellar medium. This 
list predicts a set of products which ought to be present 
in the interstellar medium, and also identifies some molecules 
which are not formed by this mechanism. We quantify this 
mechanism with a crude choice of values for the parameters 
involved, and compare the predicted abundances with those 
observed in the direction of C Oph. 

H . INTERNAL ENERGY IN CHEMICAL EXCHANGE REACTIONS 
The relative translational energy in a neutral chemical 
exchange reaction has traditionally been considered the 
important factor which may cause the reaction to proceed, i.e. 
it may assist the reactants in overcoming the activation 
energy barrier nearly always present. The significance of 
the vibrational energy in chemical exchange reactions has, 
however, been emphasized in several recent experiments. 
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For example, Jaffe and Anderson (1968) have shown that 
vibrational energy is essential to bring about a reaction 
between HI molecules. There is evidence (Parker and Pimentel 
1969, Polanyi and Tardy 1969) that in the endothermic reaction 
between H and HF the vibrational energy of the HF molecules 
is the driving energy. Other examples are given by Spicer 
and Rabinovitch (1970). Theoretical investigations by 
Polanyi and collaborators also place emphasis on the vibrational 
energy of the participant molecules. Their work suggests 
that vibrational energy is particularly effective in forcing 
otherwise endothermic reactions to proceed (Polanyi and 
Wong 1969, Mok and Polanyi 1969). 

Let us consider an example of such a reaction which 
has relevance to the interstellar medium. The enaothermicity 
in the reaction O + “* OH + H is 0.082 eV (Wilkinson 

1963, Gaydon 1968). This reaction has been suggested by 
Carroll and Salpeter (1966) as the formation mechanism for 
OH in the interstellar medium. Assuming that loss of OH 
occurs predominantly in the inverse reaction then Carroll 
and Salpeter found that in short lived high temperature 
periods - during which the rates are high - amounts of OH 
consistent with the observations were formed. However, photo- 
dissociation of OH certainly occurs, directly, and also by 
predissociations. A mean life in the typical interstellar 
field is ~ 10^s (Stecher and Williams 1966). Dalgarno and 
McCray (1972) have also pointed out the possibility of the 
fast reaction C + + OH as a loss mechanism for OH. Consequently, 
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the reaction cannot supply ..enough interstellar OH if 
translational energy is the driving energy. However, in 
the light of the discussion given above, we expect that 0 + H 2 
reaction will proceed if H 2 is vibrationally excited (written 
H 2 )with v" ^ 1, for this level is energetic enough to 
overcome the measured activation energy of 0.44 eV (Wong 
and Potter 1965, Westenberg and de Haas 1967) plus the 
endothermicity . We now assume that this is the case, and 

3f * 

that the reactions with H 2 proceeds with a rate constant 
~10 cm s" comparable with the "pre-exponential factor" 
inferred from experiments listed which are driven by 
kinetic energy. Since higher vibrational energy usually provides a 
bigger reaction rate, we are here making a conservative 
assumption. , 

Though the discussion above is concerned with O 

reacting with H„ the situation for other atoms is not very 

different. Metathetical reactions usually involve an activation 

energy rather less than 1/2 eV (e.g. Kaufman 1969). The 

endothermicity may be calculated from the binding energies 

of the molecules involved. The assumption is that if enough 

energy is inserted in the form of vibration in Hg then the 

-11 3 -1 

reaction proceeds with a rate coefficient typically ~10 cm s , 
(Kaufman 1969). 

III. THE SCHEME FOR DIATOMIC MOLECULE PRODUCTION 
We have recently estimated the abundance of vibrationally 
excited Hg in interstellar clouds (Stecher and Williams 1972b). 

The majority of absorptions in the Lyman and Werner bands 
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cascade down into highly excited vibrational levels of the 

groundstate. Because Hg is homonuclear, these levels 

6 

have long radiative lifetimes, ~10 s (Dalgarno, private 
communication), and we calculate that -0.1% of H 0 is in some 
excited vibrational level in typical clouds of H 0 . In Table 
I we show the fraction a v «i of molecules in each vibrational 
level allowing for excitation in the Lyman and Werner bands 
and vibrational relaxation in the ground state. This 
table was computed using the data of Allison and Dalgarno 
(1969) and Cartwright and Drapatz (1970). 

Since according to the discussion above, Hg is able to 
react with atoms and ions with high rates, we draw up a set 

a|c 

of allowed reactions between Hg and the most abundant atoms 
and ions in the interstellar medium. These are C* r , 0, N, Si^, 
Mg + , S + . These reactions are presented in Table II. The 
primary set of reactions is of H* with each of the atoms 
and ions. This set gives rise to hydrides and hydride ions, 
which in secondary reactions with the atoms and ions present 
produce a great variety of subsequent diatomic molecules 
and molecular ions. Further reactions between these second 
stage products are possible, but their products are expected 
to be down in abundance by a large factor. This assumption 
is later justified by numerical work. Apart from the 
primary reactions of 0 and N atoms with all other 

^ y 

reactions are of ion-molecule form and are expected to have 

very high rates (McDaniel et al. 1970). The products of 

these primary reactions may be formed in excited vibrational 

4 



levels. In many cases it is extremely likely (Kuntz et al. 

1966). However, these product molecules all possess dipole 

moments, and so vibrational relaxation will be relatively 

rapid. Reactions between neutrals, with very few exceptions, 

(Williams 1971, Kaufman 1969) involve activation energies 

which effectively prevent these reactions occurring with 

appreciable rate at the temperatures of interstellar clouds. 

We therefore ignore neutral-neutral reactions, apart 

* 

from the primary reactions with 

One of the most interesting features of the set of 
reactions in Table II is that it gives an immediate prediction 
of which molecules are produced, and which are not produced, 
by these mechanisms. The criterion is that all the reactions 
listed shall be exothermic. Endothermic reactions will have 
rate coefficients smaller by about four orders of magnitude 
in interstellar conditions. Exothermicity was ensured 
by reference to the tables of Wilkinson (1963) and Gaydon 
(1968). Reactions which do not appear in Table II have 
been omitted by virtue of their endothermicity . We list 
in Table III the obvious absentees from the products. An 
observation of any of these molecules or ions in the interstellar 
medium would immediately invoke a different mechanism. The 
allowed products can be seen from Table II. In Table III 
some of the neutral absentees could be formed by recombinations 
of electrons with suitable ions which are present. These 
molecules are noted in the table. However, if dissociative 
recombination occurs it will be a much more rapid process. 
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Typically, it has a rate coefficient ~ 10 -7 cm 3 s -1 (Bardsley 
and Biondi 1970) . The most favorable theoretical rate for 
recombination (including dielectronic recombination) is 
~ 10 -10 cm 3 s -1 , but rates of this magnitude have not yet 
been found. According to Hasted (1964) it is likely that 
recombination rates are two orders of magnitude less than 
this. We conclude that neutral molecules are unlikely to 
be formed by ion-recombination with electrons, and that any 
interstellar molecular production theory based on this 
assumption is in conflict with present laboratory measurements. 

The gratifying result of the qualitative inspection 
of these equations is that those diatomic molecules which 
have been observed in the interstellar medium are, indeed, 
formed by the processes described here. 

iv. numerical predictions of molecular abundances 

If one can assign suitable rate coefficients for the 
reactions in Table II and also for all relevant loss mechanisms, 
then one may, on the assumption of dynamical equilibrium, 
calculate equilibrium abundances. The loss mechanisms 
for hydrides, of course, involve the secondary reactions, 
but all molecules are subject to photodissociation, neutrals 
also to photoionization if their ionization potential is 
less than 13.6 eV, and ions are also subject to dissociative 
recombination with electrons. Photodissociation, photoionization, 
and the amount of vibrational excitation all depend on the 
interstellar radiation field which varies greatly from place 
to place. 
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It seems appropriate to calculate molecular abundances for 
a specific situation in the interstellar medium so that a 
comparison with observation may be made. The most careful 
and complete observations available to us at this time are 
those which were made in the direction of C Oph (Herbig 
1968, Smith and Stecher 1971, Smith 1972). We, therefore , 
attempt to define this situation so that we may compute the 
appropriate parameters. 

The star C Oph is a bright early type star (0 9.5V). 

We shall represent its spectrum in the relevent wavelength:. 


regions by that of a 30,000°K black body diluted by a 
factor 10 -14 . This implies that we are concerned with a 
region some parsecs away from the star. Herbig s detailed 
observations of C Oph coupled with Stecher' s (1968) measurement 
of Lya in this direction led Herbig to suggest that the 

2 

molecules were found in a relatively dense sheet (n(H) ~ 10 -10 



of gas in which the extinction was about 1 mag. Even if 
the situation were dynamically steady (which it is not) then 
we would expect on the basis of production theories 
considerable amounts of H 2 in clouds of this density (Hollenbach, 
Werner and Salpeter 1971). In fact, C Oph is moving at high 
velocity relative to its local environment (Blaauw 1961) 
and clouds of density ~ lO^cm ^ in a less intense radiation 
field will certainly contain much Hg. In the computations 
of photodissociation, photoionization, and n(H 2 ) we assume 
the radiation field stated above, which is equivalent to 
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~10^ photons cm - ^s~^8 ^ at 1000$ and 5 x 10^ photons 
cm”^s^8”^ at 3000$. Inside the cloud extinction reduces 
these fluxes by a factor of 5 at 1000$ and factor of 3/2 at 
3000$, using Stecher's (1969) extinction measurements. This 
field', including extinction, gives an excitation rate from 

it -9 -1 

X, v =0, excluding dissociation, of 6 = 10 s 

We now discuss the various parameters and the values 
adopted. 

(a) Reactions of H 2 with neutral atoms. 

3j( 

The case of 0 + Hg has been discussed already. We 

-11 3 -1 

adopt a rate coefficient of 10 cm s for v M ^ 1, which 
is probably a conservative estimate. No similar measurement 
involving N appears to be available, but other reactions 

between N and O with the same molecules have similar rates. 

-11 -3 -1 * 

We, therefore, adopt 10 cm s for the reaction H 2 + N, 

which becomes exothermic for v" ^ 4. 

j|C 

(b) Reactions of H 2 with ions. 

* X 

We have discussed H 2 + C elsewhere (Stecher and Williams 
1972a). The rate coefficients used were based on a 

_ Jg 2 

measured cross section of 10” cm (Maier 1966) , and the 

-11 3 -1 

rate coefficients adopted were 10” cm s . It is possible 
that this is an underestimate. The measurements were performed 
at relative energies of several volts to overcome the 
endothermicity of 0.4 eV. If this problem is removed by 
the use of internal energy, then the reactants can approach 
and react at much lower velocities, at which the usual 
polarization forces dominate (Gioumousis and Stevenson 1958). 
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One could argue that a more appropriate value would be 

one or two orders of magnitude larger. However, we retain 

-11 3 -1 

the values of 10 cm s as a conservative estimate. 

•' ■ *; . - 

+ + 

Measurements for Si and S have not been performed. We 
adopt the same value for the analogous reactions when they 
are made exothermic by the insertion of sufficient vibrational 
energy* Again, this may be an underestimate. 

(c) Other ion - molecule reactions. 

These reactions are almost invariably fast when they 

occur, i.e. when they are exothermic. Many examples are 

available (McDaniel et al. 1970, Ferguson et al. 1969) and 

-9 3-1 

nearly all have rate coef f icienls ~ 10 cm s There is 

little direct evidence concerning the actual reactions 

— 9 3 — 1 

involved here. We adopt rate coefficients of 10 cm s 
for all ion-molecule reactions, excluding those discussed 
in (b) . 

(d) Dissociative Recombination 

As mentioned earlier, this process is extremely 

—6 3—1 

rapid when it occurs. Values as great at 10” cm s“ are 

measured for large complex ions bound by polarization forces. 
—7 3—1 * 

But values ~10 cm s seem appropriate for diatomic ions. 

When neither atom is hydrogen, the complexity of states 
make it extremely probable that there is a favorable crossing 
of a neutral repulsive potential curve with the ground 
state curve of the ion. Thus, the rates will be large. 

For hydrides the situation is simpler, and the existence of 
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suitable crossings is not predictable. However, unpublished 
calculations by Krauss indicate that there is a suitable 
crossing in CH* + e so that the rate for dissociative 
recombination will be high. In the absence of further 
information we assume all molecular ions in this work 

-7 ^ 1 

dissociatively recombine with a rate coefficient of 10 cm s . 

(e) Photoionization. 

Neutral molecules will be readily photoionized if 

their ionization potentials are appreciably less than 13.6eV. 

We choose a cut-off at 13.0 eV. We assume that near 

- 17 2 

threshold the cross sections are all ~10 cm . The important 
contributors are those listed in Table IV. 

(f) Photodissociation. 

The photodissociation cross sections for the molecules 
involved are poorly known, if at all. Most effort has been 
applied to the CH molecule in recent years, and the discovery 
of many predissociated levels (Herzberg and Johns 1969) 
has led to a great reduction in the lifetime of CH against 
photodissociation in the interstellar medium (Hesser and 
Lutz 1970). Recent estimates by Stief et al. (1972) show 
that a strongly bound molecule like CO may still have a 
fairly short lifetime in the interstellar medium depending 
on the photon efficiency. None of the calculations 
are directly applicable to the present case. Correction 
must be made for the stronger radiation field we use and the mean 
extinction. We have made crude estimates of photodissociation 
taking account of electronic structure and states where 
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. _g -1 

known. We have adopted two rates ; ' fast', 10 - s * , and- slow, 
10 _9 s"- 1 -, ' as shown- in Table V. ‘ "'?• • a ■**'■> *.-• •/ 

(g) The abundances assumed were the usual cosmic v’aTaes ' f- 
(Allen 1963). We have further assumed ri (Hg/ - ;n(H) /'so 
that the abundances are calculated relative to the total 
density of n = n(H) + 2n(H 2 ) . We have adopted n(e) = 

... j.a *’ *• • ,.. t . - ^ ~’r ,v ^ ^ 

6 k 10“ n which is equal to that of n(C ), and n(H ) = 

3 x 10“ together. We expect a contribution to both e and 
H + abundances from direct photoionization of Hg for v" ^ 4 
which is permissible for radiation for which X > 9128 (Steelier 
and Williams 1972c). 


V. RESULTS AND DISCUSSION ' 

- We present the equilibrium abundances of all the product 

. -t, ' '• i, r i :r , ' . _3: ' : ■ r ' / 

molecules at a specified density of H, of 10 cm in Table VI 

■ ' ■ -,V' .-V=. - .» - ■ v 

and in Figure I we show the abundance of the' observed molecules 

as a function of n(H 2 ). Clearly, the choice ot ? the? parameters 

involved is difficult and, at times/ arbitrary. Hence, absolute 
reliance cannot be placed on. the individual abundances. Howeyer, 
the figure shows results which are quite striking’ ih thedr 

behavior. For the first time, here is a mechanisin' which 

readily produces molecules XY in greater abundance than ' 

molecules XH, where X and Y are not H. The ratios between 

observed, molecules are in harmony with the observations-, 

.••••■* ; '1 V ‘ •/? :>/ 

remarkably so, considering the crudeness' of the parameters. 

’/// ; -iv. J At .> . ; j- - - 1 -- '■ 

For example the ratio of [C0]/[CH] observed is aboiit 40 . 

The figure leads to a ratio on the order of 10 over a wide 

of densities. The measured [.CH.]/[CN] ratio is about 5. 
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For densities less than 500 cm“^ we calculate this ratio to 
be between 1 and 5. The measured [CH]/[CH + ] ratio is on 
the order of unity, consistent with our calculations. 

The absolute values of the column densities may be 

calculated using Stecher's (1968) column density for C Oph 

as a normalizing observation. Calculations are presented 

in Table VII. There is remarkable agreement between some 

theory and observations at the density of n(H 2 > = 10 cm , 

but a discrepancy for CO. Slightly higher densities 

would improve the agreement. Insufficient CO is produced, 

by more than a factor of 10. The source of this discrepancy 

may lie either in the choice of the photodissociation rate, 

or in ignorance of a particularly favorable production mechanism. 

A new measurement of CO photodissociation would be valuable. 

Within the crudity of the present work, this discrepancy is 

not too serious. The amount of OH calculated is above the 

upper limit set by Herbig. This is not a serious objection 

for several reasons: OH has yet to be seen in the UV 

though its presence is well established, there appears 

to be some observational problem; it is possible that OH is 

particularly vulnerable to third stage reactions such as 

+ OH HgO + H which should proceed with a rate coefficient 
-11 3 -1 

~ 10 cm s (Dixon- Lewis, et al. 1966). We have not 

included these third stage reactions in the computations, 
for lack of data on comparable reactions. However, the 
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reaction referred to here should constitute a significant 
loss rate for OH and an interesting source of H o 0 in the 
interstellar medium. There are some other reactions of 

. - .• j * 

a similar nature which may also be fast, but which we have 

not included in the production scheme because of lack of 

information concerning the complete set of an logous 

* 

reactions. It is possible that NH and Hg react to form NH 2 

similar to OH and Hg. Ion-molecule reactions more complicated 

than those considered here may also occur, for example 

CH + + H 2 “* CHg + H is known to be fairly fast (McDaniel et-.al, 1970) 

Thus, although we have limited our discussion to diatomic 

molecules and ions, it seems quite likely that triatomics 

will also evolve naturally. However, we need more information 

before we can include them in our computations. 

We have chosen to examine the situation in the direction 

of C Oph where we believe the radiation is intense. It 

is interesting to consider the effect of reducing the 

radiation intensity corresponding to other regions in the 

interstellar medium. As far as neutral molecules are concerned 

their loss is dominated by photodissociation, and their 

production depends directly or indirectly on the abundance 

of H 2 whioh in turn is dependent on the radiation field. 

Hence, rather surprisingly, the abundance of neutral molecules 

formed by these mechanisms is not expected to vary greatly 

with radiation field. However, the dominant loss mechanism 

for molecular ions is dissociative recombination, and the 

electron abundance as fairly insensitive to radiation intensity, 

over a wide range. Hence the abundance of molecular ions 
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should be roughly proportional to the radiation intensity. 
These simple arguments may be used as a tool to measure 
the interstellar radiation field when the rates are more 
precisely known. 
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LEGEND FOR FIGURE 

3 + 

Calculated abundances per cm for the molecules CH, CH , 

CO, CN observed in the direction of C Oph, as a function 



MOLECULAR ABUNDANCES PER 



Equilibrium fraction of H 0 in each excited vibrational level; 
in terms of the excitation rate S (s ) . In our model we set 
8 - 10 -9 s -1 


Vibrational level 


x Fraction of H 2 in level 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


0.994 
0.467 
\ 0.288 
-■ 0.198 
0 . 142 
0.105 
7 0.079 
0.061 
0.048 
0.038 
0.032 
0.028 
0.024 
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0.021 


TABLE II 


Exothermic Reactions 

Primary Reactions 
H*(v" * 1) + C + - CH + + H 
H*(v" * 9) + C + ~ CH + H 
H*(v" * 6) + Mg + - MgH + H 
H*(v" * 3) + Si + - SiH + + H 
H*(v" s 2) + S + - SH + + H 



H^v" * 1) + 0 - OH + H 
h£(v” * 4) + N - NH + H 







TABLE III 


A List of "Absentee" Molecules 



The molecules in the left hand column may be formed by 
recombination of electrons from ions which are present 
(see text) . 



TABLE IV 


-9 - 1 

Photodissociation Rates, C, in Units of 10 s for the 
Radiation Field Described in the Text, and Allowance 

for Extinction 


Molecule 


C 


CH 

2.2 

C 2 

oo 

• 

o 

SC 

2.9 

NO 

3.9 

°2 

00 

• 

o 


0.8 





TABLE V 


Photodissociation Rates Adopted for the Radiation Field 
Described in the Text, with Allowance for Extinction. 
"Fast" implies a rate of 10 - 8s~l, and "slow" a rate of 

lO-Ss" 1 . 


Fast Slow 


OH 

SN + 

CH + 

CN 

CH 

NO 

0H + 

SiO + 

NH 

no" 1 " 

CO + 

SiC + 

SiH + 

°2 

CO 

SiN + 

SH + 

°2 + 

=2 + 

MgO + 

SC + 

NH + 

C 2 


SC 

MgH + 

CN + 






TABLE VII 


A Comparison Between Observed and Calculated Column 
Densities in the Direction of C Oph. The CO 
Observation is from Smith and Stecher (1971) , the 
Other Observations from Herbig (1968). The Calculations 
Refer to the Model Described in the Text, with n(Ko)= 
lO^cm -3 . The Column Densities are in Units of cm - 2. 


Molecule Observation Calculation 


CO 

1.8 

X 

10 15 

6 

x 10 13 

CN 

8.7 

X 

10 12 


10 12 

CH + 

2.6 

X 

10 13 

3 

,,,13 
x 10 

CH 

4.3 

X 

10 13 

0.8 

1(1 i3 
x 10 

NH 

<7.4 

X 

10 12 

8 

in 12 
x 10 

OH 

<8.3 

X 

io 13 

2 

m 14 

x 10 

MgH 

<4.3 

X 

10 12 


0 

SiH 

<5.6 

X 

10 12 


0 

CO + 

<3.5 

X 

10 13 

3 

x 10 11 


NASA— GSFC 


